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Abstract 

We generalize the concept of Fermi normal coordinates adapted to a geodesic to the case where 
the tangent space to the manifold at the base point is decomposed into a direct product of an 
arbitrary number of subspaces, so that we follow several geodesies in turn to find the point with 
given coordinates. We compute the connection and the metric as integrals of the Riemann tensor. 
In the case of one subspace (Riemann normal coordinates) or two subspaces, we recover some 
results previously found by Nesterov, using somewhat different techniques. 
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I. INTRODUCTION 



The construction of Riemann normal coordinates is well known. For any point p of a 
Riemannian or Lorentzian manifold {Ai,g) and any vector Vp at p there exists a maximal 
geodesic 7v(A) with starting point p and initial direction Vp. We define the exponential 
map expp that takes a subset of Tp, the tangent space to at p, into Ai, such that expp(V) 
is the point q a unit parameter distance along the geodesic 7v from p. 

We can choose an orthonormal tetrad basis {E(a)} at p and then define the coordinates at 
q by the relation q = exp(x°E(Q,)). Such coordinates are called Riemann normal coordinates. 

The Fermi normal coordinate construction [l| is also well known. We start with a timelike 
geodesic 7k (A) with tangent vector K at p. (We will consider only geodesies, not arbitrary 
timelike curves.) Given any vector V G Tp{Ai), we can write V = A + B, where A is in 
the direction of K and B perpendicular to K. We then let q = expp{A) and define a map 
Fermip such that Fermip(V) = expq{B). That is to say, Fermip(V^) is found by first moving 
unit distance along the geodesic 7a from p to q. We parallel transport B from p to q and 
then move unit distance along the geodesic whose tangent vector at g is S. 

We can define an orthonormal basis {E(q,)} at p such that E(o) is parallel to K. Then 
A = x''E(o), B = x*E(j), giving the usual construction of Fermi normal coordinates 

In this paper, we will generalize this construction to allow an arbitrary number of arbitrary 
subspaces, rather than just a timelike geodesic and the perpendicular space, and an arbitrary 
number d of dimensions. In Sec. [TTl we will construct the generalized coordinate system, in 
Sec, mil we compute the connection, and in Sec. II VI we compute the metric in the generalized 
Fermi coordinates. We conclude in Sec. El 

We use the sign convention (+, +, +) in the classification of Misner, Thorne and Wheeler 

0. 

II. MULTI-STEP FERMI COORDINATES 

Consider a ci- dimensional Riemannian or Lorentzian manifold {Ai,g). We will start 
our construction by choosing a base point p ^ Ai. We decompose the tangent space Tp 
into n subspaces, Tp = A^^ x A'^p^ x A'^^ . . . x A^^ so that any V G Tp can be uniquely 
written as V = V(i) + V(2) + V(3) + . . . + V(„). We choose, as a basis for Tp, d linearly 
independent vectors {E(q,)} adapted to the decomposition of Tp so that for each m = 1 . . .n, 

{E(q) |a G Cm} is a basis for A^\ where ci, C2, . . . c„ is an ordered partition of {1 . . . d}. Thus 
each V(m) = J2a€cm ^"^-^(a)- ^hc vectors {E(a)} need not be normalized or orthogonal. 

The point corresponding to coordinates x° is then found by starting from p and going 
along the geodesic whose whose tangent vector is V(i), parallel transporting the rest of the 
vectors, then along the geodesic whose tangent vector is V(2), and so on. An example is 
shown in Fig. [TJ With that general construction of multi-step Fermi coordinates we can 
define a general Fermi mapping q = Fermip (V) given by 

9(0) = P 

9(1) = expp(V(i)) 

9(2) = expg^^^(V(2)) (1) 

1 = Q{n) = exp (V(„)) 
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FIG. 1: Construction of 3-step Fermi coordinates in 3 dimensions. We travel first in the direction 
of E(]^), then E(2), then £(3-), parallel transporting the triad as we go. 

From that general construction we can return to the original Fermi case by choosing ci = {t} 
and C2 = {x, y, z}. In the Lorentzian case, we could also choose a pseudo-orthonormal tetrad 
E„, E^, Ej;, Ey, with and E^ null, E^ ■ E„ = —1, and other inner products vanishing, and 
Co = {u} and ci = {v,x,y}. 
For later use we will define 

m 

V(<^) = E ^"E(.) = E^(0 (2) 

Q£ClUC2U---UCm 1 = 1 

^{<m) = ^(<(m-l)) = ^{<m)) " ^ (m) (3) 

Then we can write = Fermip(V(<m,)). 

An example of a spacetime where these multi-step Fermi coordinates might be used comes 
from brane-world models. A general brane-world metric with one extra dimension is 

ds^ = b{w) [-dt^ + a{t) {dx^ + dy"^ + dz^)] + dw"^. (4) 

A simpler metric of that form is used, for example, in j3|. In this kind of spacetime it 
might be useful to introduce three-step Fermi coordinates with: Ci = {w}, C2 = {t} and 
C3 = {x,y,z}. 



III. CONNECTION 

We will parallel transport our orthonormal basis vectors E(„) along the geodesies that 
generate the coordinates, and use them as a basis for vectors and tensors throughout the 
region of A4 covered by our coordinates. Components in this basis will be denoted by Greek 
indices. We will use Latin letters from the beginning of the alphabet to denote indices in 
the Fermi coordinate basis. Of course at p, there is no difference between these bases. 

Latin letters from the middle of the alphabet will denote the subspaces of Tp or equiva- 
lently the steps of the Fermi mapping process. 

We would like to calculate the covariant derivatives of the basis vectors, V^E(a). Having 
done so, we can calculate the covariant derivative of any vector field V = l^'^E(^) along a 
curve /(A) as 



dX dX \dX " 
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FIG. 2: The covariant derivative V^E(q,) is the change in E^^) under parallel transport along the 
path qi ^ p ^ q2 ^ qi, divided by dx. The parallel transport can be decomposed into a series of 
transports clockwise around trapezoidal regions with sides AE^^^dx and XdA. 

To evaluate V/3E(q,) at some point qi = expp(X), consider an infinitesimally separated 
point q2 = expp(X + E(^)(ia;). The covariant derivative of E(q,) at qi is the difference between 
E(Q,)(g2) parallel transported to qi and the actual E(Q,)(gi), divided by dx. That difference is 
the same as the change in by parallel transport around a loop following the geodesies 
from qi backward to p, the infinitesimally different geodesies forward from p to q2, and the 
infinitesimal distance back to qi. We can write this loop parallel transport as an integral 
over the Riemann tensor. 

Let us first consider the Riemannian shown in Fig. [2J The total parallel transport 

can be written as the sum of parallel transport around a succession of small trapezoidal 
regions whose sides are AE(/3) and XrfA. By using the definition of the Riemann tensor we 
have ^ 

VpEj.= [ dXR\si3{XVXX'. (6) 
Jo 

Here R is evaluated at the point expp(AX), which we have denoted merely AX for compact- 
ness. 

Equation reproduces Eq. (13) of Ref. [3]. Note, however, that Eq. (jS]) is exact and 
does not require R to be smooth, whereas that of Ref. 0] was given as first order in R and 
was derived by means of a Taylor series. 

We see immediately that the covariant derivative of any E(q,) at X in the direction of X 
vanishes. This happens simply because changes with dX. in the direction of X correspond 
to additional parallel transport of E(q,). 

Let us now consider the general case where there are n steps, and compute V^E(q,). Since 
the coordinates are adapted to our construction, the index /3 must be in some specific set 
Cm, which is to say that the direction of the covariant derivative, E(^), is part of step m in 
the Fermi coordinate process. We will write the function that gives that m as m(/3). Some 
particular cases are shown in Fig. |3l 

If m = n (leftmost in Fig. |3]), only the last step is modified. The integration is exactly 
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m=n=3 





P •' P P 

FIG. 3: Original and displaced geodesies for Fermi coordinates with n = 3 and m = 3, 2, and 1. 

E(p)dx 




FIG. 4: Part of the calculation of V pE(^^-^ in the case n = 2, m(/3) = 1. The geodesic of the 
first step has been modified, causing the geodesic in the second step to be displaced. The parallel 
transport integrates the Riemann tensor over a series of rectangular regions between the 2 second- 
step geodesies. 

as shown in Fig. O except that it covers only the final geodesic from X(<„) to X(„), 

= / ^'^-R'^a<5/3(X(<„) + AX(„,))AX(^„). (7) 

If m < ra, then we are modifying some intermediate step, and the path followed at later 
steps is displaced parallel to itself. In that case we get an integral over rectangular rather 
than trapezoidal regions, as shown in Fig. |H For general m there is a contribution for each 
step j > m. The j = m contribution integrates over trapezoids that grow with A, while 
the j > m contributions integrate over rectangles with fixed width dx. We can write the 
complete result 



n „i 

where m = m{f3) and 

'^jm(A) — 




(8) 



(9) 
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Equation ([8]) is exact and includes Eqs. ( 16|7|) as special cases. 

Consider the case where ci consists only of one index. If m > 1, there is no j = 1 term in 
Eq. ([8]). If m = 1, then (3 is the single index in ci, and X^^-^^ = unless 5 = /3, so the j = 1 
term in vanishes because Ra-yS/s is antisymmetric under 6 ^ (3. Thus there is never a j = 1 
contribution to Eq. ([S]) when there is only one index in Ci. 

Now suppose X lies on the first generating geodesic, so X(j) = for j > 1. Then all j > I 
terms vanish in Eq. ([H]). So if Ci consists only of one index, all Christoffel symbols vanish at 
X. This is well known in the case of the usual Fermi coordinates. 



IV. METRIC 

Now we would like to compute the metric g at some point X. Specifically, we would like 
to compute the metric component Qab in our generalized Fermi coordinates. 

We will start by considering the vectors Z(a) = d/dx"". These are the basis vectors of the 
Fermi coordinate basis for the tangent space, so the metric is given by Qab = T^{a) ■ Z(b). Thus 
if we compute the orthonormal basis components Z^"^ we can write gah = Va/sZ^^-^Z^j^y 

Again we will start with the case of Riemann normal coordinates. Let W{t, s) be the 
point exppS(X + tE(„)). Define Y = dW/dt and V = dW/ds. Then Y(X) = Z(„) and 
y/^ = X^ + t6^. The components of Z(„) at X can be calculated by integration, 

zf„,(x) = r^(x) = ^'d.^^. (10) 

Because the orthonormal basis is parallel transported we can write 

By construction, the Lie derivative LvY = and thus Ch. 4] 

DY _ DV 

ds dt 

From Eq. ^ we have 

DV^ dV^ 



(12) 



derivative is with respect to the orthonormal basis. 
From Eq. ([6]) we have 



+ l^^F"V«Ef^^ = 5f + sd^V^W^El^ + 0{R^). (13) 

(7) 



dt dt ^ " W " ^ " " W 

where we have retained 6^ instead of writing Va-E'fL^ to make it clear that the covariant 



VaE^JsX)= [ dX\R'^^sa{\sX)sX^ = - [ c/AAi?%„(AX)X'^ (14) 
Jo ^ Jo 

Taking t = 0, V is just X. Combining Eqs. (I10til4p . we find 



^a+S2 [ rfAAi?%,(AX)X^X^ 
Jo 



+ 0{R^ 



+ ! d\\{l-\)R^^sa{XK)X^X^ + 0{R^). (15) 
Jo 
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From Eq. ( !T5|) . the metric is given by 

gab = Vab + 26:6^, f dX A(l - X)R^,sp{XK)X'X^ + 0{R^). (16) 
Jo 

Equation f|T6|) reproduces Eq. (14) of Ref. 0]^. 



Next let us consider the case where there are n steps in our procedure. We will define a 
set of functions Wj as 

Wj{s) = Fermip(X«,) + sX(,)). (17) 
The path Wj{s),j = 1 . . .n, s = . . .1 traces the geodesies generating the Fermi coordinates 
for the point X. Now consider Z(a) = d/dx"'. Let m = m{a), so Z(a)(p) G A^K Then let 

{sX(j) j <m 

X«j) + s(X(j) + tE^a)) J = m (18) 
X«j) + tE(„) + sX(j) j > m 

Let Y = dW/dt and V = dW/ds as before. To find Z(a) we now must integrate over a 
multi-step path from p, 

" 9y''(WG(s)) 



ds 



^(')W = E/ 

The generalized version of Eq. f lT3|) is 

j < m 



(19) 



Dl^/'(W^,(., t)) _ rfy/' ^ y7y«v„<, = 5f + .5."l-W.<, + 0{R^) J = m (20) 



dt dt ^ " W l "^T._ i^,^ 

Now 



VaE^^^{W,{s)) = J2 ^TT / dXakm{X)R^M^«k) + AX(,))Xf,) (21) 



fc^^^fcil^J Jo 

where 



..,(.) = <!j (22) 



The k = j term is analogous to Eq. (jT4l) . while the others have no dependence on s. 
Combining Eqs. (HHnUHHSI]) we get 



^(a)(X)=5f + Fi^ + 0(i?^) (23) 



where 



F^ = J2J2 ' rfAafc„(A)i?^,,(X«,) + AX(fc))X(i)Xj.) 

n i /.I 

= Y.Hj d\akrn{X)hkj{X)R^M^{<k) + >^'y-{k))Xl)Xl^ (24) 

^ Ref. uses the same sign convention for R^fs^s as the present paper, but the opposite convention for gat 
and consequently also for Ra/s-yS- 
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where m = m{a) and 

6fc,(A) = <;: ^ (25) 

So the metric is 
where 



gab = VafiZ^a) 4) = ^^^^ + + F,a + 0{R^) (26) 



n 3 pi 



(27) 



j=m k=m 



where m = m(/3). 

Once again consider the case where Ci contains only one index. As discussed with respect 
to Eq. dHD, if /3 e Ci, there is no nonvanishing k — 1 term in Eq. 0271) . Thus Qat — f]ab ci-t 
points on the first generating geodesic. This is also well known in the usual Fermi case. 

Now suppose Ci consists only of one index and furthermore n = 2. The only possible 
term in Eq. fl271) is then j = k = 2, so 



FaP = / '^Aa2m(A)(l — A)-Ra^5_g(X(i) + AX(2) )X(^2)-^(2) ■ 




(28) 



where m = m (/3). Equation Eq. ([28D is equivalent to Eq. (28) in Ref. 8 in the case where 
the generating curve of the Fermi coordinates is a geodesic. 

Now we are in a position to discuss the region of the manifold over which the multi- 
step Fermi coordinates are well defined. Assuming there are no singularities or edges in 
the manifold, the only problem would be if the metric Qab is degenerate, which in turn can 
happen only if Fab is of order 1. Thus the Fermi coordinates will be well defined providing 
that (3 

\Ra-y^5\{X'f -^1 (29) 

throughout the region of interest, for all a, 7, (3, 6, e. 

In the case where there is only one index in ci, there is no contribution to Fab from 
X(i). Then it is sufficient for Eq. (129!) to hold for e > 1. In other words, if the first step is 
one-dimensional, it can be arbitrarily long pL]- 



V. CONCLUSION 

We have generalized the usual Fermi normal coordinates in the case where the generating 
curve is a geodesic to allow for any number of steps and for a subspace of any dimension 
at each step. We have derived the connection (exactly) and the metric (to first order in 
the curvature) as integrals over the Riemann tensor. Our results reproduce several formulas 
previously derived by Nesterov 8 with a more geometric approach in a more general setting, 
and without reference to any derivative of R. 
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